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I. INTRODUCTION 
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We study slow-roll inflation with the Gauss-Bonnet and Chern-Simons corrections. We obtain 
general formulas for the observables: spectral indices, tensor-to-scalar ratio and circular polariza- 
tion of gravitational waves. The Gauss-Bonnet term violates the consistency relation r = — 8nT. 
Particularly, blue spectrum wt > and scale invariant spectrum ISnTl/T" <C 1 of tensor modes 
are possible. These cases require the Gauss-Bonnet coupling function of ~ 10*/Mpi. We use 
examples to show new-inflation-type potential with lOA/pi symmetry breaking scale and potential 
with flat region in > lOMpi lead to observationally consistent blue and scale invariant spectra, 
. respectively. Hence, these interesting cases can actually be realized. The Chern-Simons term pro- 

^SJ ' duce circularly polarized tensor modes. We show an observation of these signals supports existence 

. of the Chern-Simons coupling function of w,^ ~ 10*/Mpi. Thus, with future observations, we can 

tr^' fi^ or constrain the value of these coupling functions, at the CMB scale. 

■ PACS numbers: 

o 

u 

• , Currently, the concept of slow-roll inflation is widely accepted by cosmologists, due to many precise observations, 
' such as the Wilkinson Microwave Anisotropy Probe (WMAP)[l[. We can explain all observations simply, using single 
II field slow-roll inflation. However, future observations might detect non-standard signals, which will be evidence for 
O a non-standard inflationary theory. The energy scale of inflation is considered to be very high compared to current 
^ \ experimental scale. Thus, we can say that a new inflationary theory might contain information with regard to high 
^ ■ energy physics, such as a superstring theory. Therefore it is important to study theories beyond standard single field 
I I ■ slow- roll inflation. 

One of candidates for a modified inflationary model is inflation with the Gauss-Bonnet and Chern-Simons correc- 
' tions. These are the only two meaningful combinations of second order curvature terms, as a low energy effective 

^ ', theory [s^l, and can also be derived from some superstring models [ll. Cosmology with the Gauss-Bonnet term can 
' affect background evolution and result in a non-singular cosmological model0-[l^ or other interestin g so lutions p!6l- 

\ [2^. In infiation with the Chern-Simons term, primordial gravitational waves are circularly polarized |26l - [30| . which 

■ is considered as an interesting target for some future observations [Sll - lssj . Many authors studied Gauss-Bonnet and 
Chern-Simons modified inflation. However, slow-roll inflation with the Gauss-Bonnet terms is not much studied, 
except in the papers [ssl [36| . Taking current observations into consideration, slow-roll infiation is the most plau- 

' sible paradigm for inflation. Therefore investigating slow-roll infiation with the Gauss-Bonnet and Chern-Simons 
^-H ,. corrections is an important subject. 

IV ' We have studied this topic in the previous paper (sBj. In this paper, we re-formalize slow-roll infiation in gravity 
. ^ ] with the Gauss-Bonnet and Chern-Simons corrections. We derive formulas for the observables, namely the scalar 
' spectral index ns, the tensor spectral index nx, the tensor-to-scalar ratio r and the circular polarization ratio H, in 
?H ] our inflationary model. The Gauss-Bonnet term violates the consistency relation r — —Sut- We show if this violation 
. . . i is observationally confirmed, the derivative of the Gauss-Bonnet coupling function at the observation scale, namely 
the Cosmic Microwave Background (CMB) scale, is fixed. We study the two typical cases, blue and scale invariant 
spectra of tensor modes. Because blue and scale invariant mean nx > and jSnxl/r <C 1, respectively, these cases 
strongly violate this consistency relation. If these are confirmed in future observations, the Gauss-Bonnet coupling 
function must be ~ 10^/Mpi, at the CMB scale. Of course, even if no such effects are observed in future, we can 
get the constraint ^ ^ < 10^/Mpi, at least. However, it is unclear whether these interesting cases are consistent with 
current observations. Therefore, we perform calculations in some concrete examples. We show new-inflation-type 
potential with symmetry breaking scale of lOil/pi leads to an observationally consistent blue spectrum; potential 
with flat region in (j) > lOMpi leads to a consistent scale invariant spectrum. In both cases, the almost constant 
is required. The Chern-Simons term leads to circular polarization of gravitational waves. We show if this circular 
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polarization is detected, we can fix the derivative of tfie Cfiern-Simons coupling function, at the CMB scale. It might 
be w ,^ ~ 10*/Mpi. Therefore, the value of the Gauss-Bonnet and Chern-Simons coupling functions at the CMB scale 
will be fixed, or constrained, with future observations. 

We organize this paper as follows: In section II, the action we considered is shown and slow-roll inflation in this 
action is studied, calculations of perturbations in slow-roll inflation is in section III, interesting features of our model 
are mentioned in section IV, and section V is for conclusion. 



II. SLOW-ROLL INFLATION 



We consider the following action: 



1 



^ J d^xV^mRoB + d''x^u:{4>)RR , (1) 



where Mpj = I/SttG is the reduced Planck mass, and RR are the Gauss-Bonnet and Chern-Simons combinations. 
These combinations are defined by 

-^GB = Ra^-fSR"'^''^ — '^RapR"'^ + R^ (2) 

RR = ^^e'^f'^'R^pp^Ry , (3) 

where e""!^^^ is the Levi-Civita anti-symmetric tensor. Hence, the Chern-Simons correction represents parity violation 
in gravity. This action is composed of the Einstein-Hilbert term and a canonical scalar field which couples to the 
Gauss-Bonnet and Chern-Simons terms through the coupling functions and oj{4i). Note that only these two 
types of higher curvature corrections, namely the Gauss-Bonnet and Chern-Simons terms, are essential as a low 
energy effective theory ^34j. Here, we take these terms as small corrections. In principle, coupling functions must 
be fixed by quantum gravity, such as a superstring theory. However, it is impossible to determine these coupling 
functions, because we haven't succeeded to construct quantum gravity. Therefore, we treat these as free functions 
and attempt to fix or constrain these coupling functions from cosmological observations. Through these attempts, we 
may get new information with regard to quantum gravity, because these coupling functions are related to properties 
of quantum gravity. 

Here, as a background metric, we choose the flat Friedmann-Robertson-Walker(FRW) metric: 

d^s = -df + a^{t)6,jdx'dx^ , (4) 

where a{t) denotes a scale factor which describes cosmic expansion. Variating Eq.([T]) in the FRW metric, Eq.(|4]), we 
get background equations: 

iMl.H^ = l<j>^ + V + ^H^i (5) 

4> + 3H,p+'^H^{H + H^)i^ + V^=0 , (6) 

where a dot denotes a time derivative and the Hubble parameter H is defined by iJ = a/ a. Because of parity symmetry 
of Eq.dH), the Chern-Simons correction cannot affect these background equations. We take < throughout this 
paper. 

Here we focus on the case in which scalar field is slowly rolling and its friction term is dominating Eq.®. In other 
words, we take slow-roll approximations: 

< 1 , < 1 . (7) 



In this paper, we consider the Gauss-Bonnet term as a small correction. Hence, it must be much smaller than the 
Einstein-Hilbert term, energetically speaking. The following inequality, hence, must be satisfied: 

H^i Hi 



2M2j7J2 2A/pj 
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Under Eqs.([7]) and ([8]), the background equation, Eq.(l5]), become 



This equation doesn't contain the Gauss-Bonnet correction since it is energetically negligible, of course. The time 
derivative of the Hubble parameter becomes 



The another background equation, Eq.([6]), is 



(10) 



(11) 



The last term in the right-hand side can be written as 

^ QUA ^ 

Hence, this term is negligible compared with the first term in the right-hand side of Eq.®, and Eq.® becomes 



= _Mp,-^ - ^ . (13) 



MpiH V 6Af3 



PI 



This equation shows that the Gauss-Bonnet coupling function ^ works as effective potential. In our inflationary model, 
slow-roll equations are Eq.([9]) and (fT3|) . We can see only the equation of motion, Eq. (fT3| . is modified. Here, we define 
the new function ^: 



C(^)-K. + S^- (14) 



'PI 

Eq. ([T3l) becomes 



MpiH V 

In our model, the expression for the e-folding number N is written as 

° , M 1 /"^ , ,y 1 



C • (15) 



N = / dtH = / d6^ = / d^— = / d^ ■ ■ -r , (16) 

where we assume that potential has a minimum at (j) — 0. 

We must consider consistency conditions for Eqs.Q and ([T5|) . First, we check slow-roll conditions, Eqs.Q: 

i2 /■2 

- M^,^ « 1 , (17) 



Second, the condition for smallness of the Gauss-Bonnet term, Eq.®, is 

?^* = -o-frCC0 = -77T7rC^7^(C-^<0)=A4i^-^4iT72 ^1 ■ (19) 



Here, we define three slow-roll parameters: 



_M|iC^0 ,_,.20 
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It is clear that if |a|,/3, 17I <C 1 are realized, all consistency conditions are satisfied and slow-roll inflation occurs. In 
the case of conventional slow-roll inflation, ^ — is realized. Hence, a and /3 becomes conventional e, and 7 becomes 
77. For later convenience, we derive expressions for H/H^, 4> and iJ^/Mpj in slow-roll parameters. From Egs.dTUl). 
and HH), we get 

H_ IV, M^,CV,_ ^ _ Hi _ 

Remember that we are taking < 0. As in a conventional scenario, the time derivative of the Hubble parameter is 
small. In the case that a can be treated as a constant, we can solve the equation for H easily: 

H = -. a = ao\t\'^^, (22) 
at 

where qq is a constant of integration. Because we choose 77 > 0, i > for a > and t < for a < 0. 

The most important difference between our model and a conventional inflation is that, in our model super inflation 
can be realized, while in a conventional one it cannot. In our model, the slow-roll parameter a can take a negative 
value. From Egs. d^ . we can see this leads to super inflation, H > 0. We can expect that if this super inflation phase 
corresponds to the CMB scale, an observationally interesting signal could be detected in future experiments, and this 
signal might be regarded as evidence for existence of gravitational higher curvature terms. 

Here, we show an interesting effect of the Gauss-Bonnet correction in slow-roll inflation. Let us consider general 
potential V {(/)), which does not satisfies conventional slow- roll conditions. If we set 

- ^ • <^^' 

C((/)) = is realized. This makes all slow-roll parameters zero, and inflation occurs. This shows that, even if we use 
extremely steep potential, slow-roll inflation can be achieved due to the Gauss-Bonnet correction. 



III. PERTURBATIONS 



In this section, we perform perturbative calculations on the background of slow- roll inflation. First, we study scalar 
and tensor modes of perturbations. Note that no vector modes plays an important role since there is no source 
for vectors. Second, we derive expressions for reconstructing the value of potential and coupling functions using 
observables. Finally, we briefly comment on the Lyth bound in our model. 



A. Scalar perturbations 



Here, we consider scalar perturbations in the flat FRW background. The metric is 

ds2 = a2(T) [-(1 + 2A)dT^ + B\,dTdx' + ((5,y -I- 2iP5,j + 2E\,j)dx'dx^] , (24) 

where r is the conformal time dr = dt/a and a bar denotes a spacial derivative. Latin indices run from 1 to 3. In the 
case of slow-roll inflation, r becomes 

^ ^ -i|tr'/" = ^-V- = ^ (25) 



aol — 1/a 1 — 1/aa 1 — a aH ' 

where r < 0. We fix gauge as E = 6(j) — 0, where Scfi is a perturbation of scalar field (j). This gauge is useful in 
this case, because the perturbative quantity ?A can be treated as a master variable and —'ip coincide with the gauge 
invariant comoving curvature perturbation TZ. Of course, this is a complete choice of gauge. Here, we concentrate on 
^p, because what we need is the power spectrum of -ip in super-horizon scale. Expanding tp in Fourier modes. 



the action for ijj becomes 



^=11"^^ I (Sf^' " Clk'M , (27) 
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where a prime denotes a conformal time derivative and 



1 - a 12 



1 - 3cr/4 ) \ M^^m 8 1 - cr/2 



H 



16 16Af2j 2i72 



(T 



and (7 is defined by 



(28) 
(29) 

(30) 



This action contains only ^. This is because scalar perturbations in the flat FRW background have parity symmetry. 
Let us define the new variable ^'k, by 4'k = A^ip\^, and the above action becomes 



S='-Ur 



This leads to the equation of motion; 



A. 







(31) 



(32) 



Here, we impose slow-roll approximation, and derive the expressions for and up to first order in slow-roll 
parameters. Remember that the time derivative of the Hubble parameter H is first order. In computing primordial 
power spectra, we treat slow-roll parameters as constants of time. We get the following equations: 



cr = 4(a - P) , 



= 2/3 



& H 



(7 = 



and 



— 1 



Al = (1 + a/2) 2/3 = 20^/3 
Let us consider the effective potential A'^/A^. For this purpose, we derive some equations: 



2C 2V 2Q 



2V, 



HP HC HV C H V H 



= -2M, 



pr 



V 



2M, 



PI" 



y2 



4a - 27 , 



A^ 



2 Al 



aB , ^ 11 + 2a - 7 1 + 3q!-7 

ai/ + = aH{\ + 2a - 7 = ^ = ^ 

I p — T 1 — a — T 



Hence 



A" 

A^ 



where is defined by 



dr 



4' 

A^ 



1 + 3a - 7 1 + 6a - 27 _ 2 + 9a - 37 _ i^V' ^ 1/4 



-f 3a — 7 



(33) 
(34) 

(35) 
(36) 

(37) 
(38) 



Eq.(|32j) becomes 



1/4 







(39) 



To quantize the quantity ^I', we replace ^I^ with the operator ^E* and expand this as 

k = Wk(T)dk + i'lk(T)aLk > 



(40) 
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where Ok and a^j^. are annihilation and creation operators. The mode function Vit.{T) obeys the same equation as ^'k 
does, namely Eq. ([5^ : 

<+(^fc'-^^^^^«k = 0. (41) 

For solving this equation, we must impose an initial condition. Because effective potential vanishes in asymptotic 
past, we can choose the Bunch-Davies vacuum at r — > — oo: 

„.^-i=e-^ (42, 

Under this condition, Eq. (PT|) has the analytic solution 



vUr) = ^e''('^/4+-'"^/2)HW(-fcr) , (43) 



where H^J^ is the first kind of Hankel function. 



= I d(logfc)^^fc^^ = / d(logfc)7's(fc) . (44) 



Here, we derive the expression for the power spectrum of scalar perturbation Vs in super-horizon, which is defined 

by 

In super- horizon, —kr <^ 1, v\^{t) has the asymptotic form 



Hence Vs becomes 



Vsik) = flMV fz^V"" = ^ f— V f^V"" (46) 

'^^"^ Ml,,A^^AlT^\v{i/2)) \ 2 ; \ 2 ) ' ^^^^ 

The scalar power spectrum is of order iJ^/ (/JMpj), and this is easily understood because Vs ^ (iJ^/0)^ ^ iJ^/(/3Mpj), 
in which we use Ea.(|2ip. The scalar spectral index ng becomes 

ns-l = 'i-2v^ = -6a + 27 . (47) 

If we take ^ = 0, the slow-roll parameters a and 7 coincide with the conventional ones e and 77, respectively. Hence, 
in this case, the expression for ng takes ordinary form. 



B. Tensor perturbations 

In this subsection, we compute perturbations of tensor modes. The metric is 

ds^ = a^(T) [-dr^ + {5,j + /iy)dxMa;^] , (48) 

where hij denotes a transverse-traceless tensor on a constant-time hypersurface. Hence it satisfies hi" — h^j]^^ = 0. 
We expand hij , which is gauge- invariant because tensor modes have no gauge freedom, in Fourier modes and circular 
polarization: 

where Pk ij ^^e polarization tensor for circular polarization and ± denotes the helicity of each polarization mode. 
Note that \/2 and Mpi are for later convenience. The action for 0^ becomes 

^ = E ^ / (^^T my? ~ c^mn , (50) 
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cr 1 k „\ „o a (J £ a , , 



where 

"""" ^ 2^w^vj ' ""^""^ ^2:2 " 2M2jA2: ■ 

The mass scale Mc corresponds to the UV cut-off scale of our model and il is defined by 

n^l^^. (52) 

2 Mpi Mpi ^ ^ 

Note that we omit ± symbols in the expressions for and for simplicity. The Chern-Simons coupling function lj 
doesn't appear in the background and scalar perturbations, but, at the same time, this appears in tensor perturbations, 
because circular modes clearly violate parity symmetry. 

From Eas. (|?T|) . we can see that if \n\ > 1, takes a negative value in one of helicity modes at the cut-off scale 
k/a = Mc- This is very important because it implies the appearance of ghost modes in gravitational waves. Since 
we don't know proper treatment for ghost modes, we impose the condition < 1. For simplicity of following 
calculation, we take as a constant, which is justified in many models because is slowly rolling. We can always 
treat the Chern-Simons correction kVl/Mca as a small correction, because |r2| < 1 and k/a < Mc- Of course, |ri| > 1 
model might lead to very drastic results and such model is interesting itself. However, we don't analyze this situation 
in this paper. 

Introducing the new variable /ik = ^t^^ > Eq. (|50p can be rewritten as 



± 

Hence the equation of motion for /j,k is 



^ = E ^ /d-/ (imLP - C^Tk'W? + . (53) 



~ 4^ ) /^ki = , (54) 



A 



where we also omit ± in a /ik expression. 

Considering slow- roll inflation, we compute spectra up to first order in slow-roll parameters and J7. The quantities 
At and Ct become 

A| = aM l-2a + 2/3T-^ j , Ct = 1 + 2a - 2/3 , (55) 



A, _ . (A|) _a(^^^ k^\ ^ ^ ^l±^±kJLn. (56) 



and we get 

At ~ 2 A^ ~ 2 V"" ^ Mctt^^ J ^ ^ "'2Mc" ^ -t ^ "'2M, 
Effective potential becomes 

_ d fA^\ fA^Y_l + a l + 2a k H 2 + 3a k H v^-l/A k H 
M-d^\A^)^\M) ^^^^^^ ~Wc-^^^~Wj^- t2 ^~Me ' ^ ' 

where is defined by 



Hence, Ea. ([M| becomes 



i^T — - + a . (58) 



^'^ + ( c'^k' ^ A ) = . (59) 



Let us quantize perturbations. We promote to the operator /ik and expand in the annihilation and creation 
operators and aly.: 



Ak = Mk(T)dk + wlk(T)aLk 



(60) 
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The equation for Uk(T) is the same as that for ^-kir), Ea. (l59l) : 



u'^+{c^^e.±J^^±^n]uk = 0. (61) 



As in the case of scalar modes, effective potential vanishes at asymptotic past. Thus we can use the Bunch-Davies 
vacuum as the initial condition for Eq. (pT|) . We take at r — >■ —oo: 

Hence, the solution for Ea. (|5T|) is 

yk = e-'^-'=-(~2CTfcr)''-y=7e-'(-/4+--W2)uQ + j.T^_l_Af7, 1 + 2i.t, 2zCTfcT^ exp (^±^^f7^ , (63) 
where U is the confluent hypergeometric function. This solution has the asymptotic form in super-horizon scale: 

"•^^^^^ V2(-CtA:t)3" r(3/2)l, 2 ) """P 1,^ 4Ct Mc j ^ ^ 

We define the power spectra for each polarization mode and the total power spectrum Vt- 

{0\hi;h^^\0) =Y,j d(logfc)^-|fc3 M! ^Y^j d(iogfc)7^± = j d(logfc)7'T . (65) 



In super-horizon, these and Vt: become 



T - M^, C| ^2 ^1^2 r(3/2) J { 2 J P 2Ct Me 



, 1 ± 17 (66) 

p,.^,l-3o + «(-^)'(^)"\ (67) 

The tensor spectral index nx is 

riT = 3 - 2i/T = -2a . (68) 

This expression becomes the same as an ordinary one in the ^ — >• limit. Remember that in our model, a can takes 
a negative value in which super infiation, _ff > 0, is realized. Hence, this expression shows that in super infiation, 
a blue spectrum for tensor modes, nx > 0, is realized. The detection of this blue signal might be evidence for this 
inflationary model, because ut < is always satisfied in a conventional model. 

We have derived both scalar and tensor spectra. Thus, we can get the observationally important quantity, namely 
the tensor-to-scalar ratio r: 

r = ^ = 16/3 . (69) 

Note that /3 — >■ e in the ^ — >■ limit. There is the consistency relation r = — Stt-t in ordinary single field slow-roll 
inflation, or r < — Snx in multi-filed one. However, these consistency relation doesn't hold in our case, because a and 
/3 are in principle different quantities. The violation of this relation might be an observationally interesting feature of 
our model. Note that a blue spectrum of tensor modes can also be interpreted as a strong violation of this relation. 

The another important feature of our model is a violation of parity symmetry. Eq. (|66p shows that the parity 
violating Chern-Simons correction results in difference between power spectra of each helicity mode. We quantify this 
difference using the ratio H which is defined by 

H= -^T ^ ^JLn . (70) 

2 Me 
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This quantity 11 takes naturally a small value, because f2 < 1 and H < Mc- However, we can still expect 11 is of 
order a few percent, which is considered to be observable in future experiments [sTI-issj . Note that we have got no 
observational constraint on ui and fl, because these only appear in tensor modes, which we haven't detected. 

In a conventional model, the transformation of potential V — >■ pV, where p is constant, has no effect on spectral 
indices, tensor-to-scalar ratio and e-folding number. Hence, an overall factor of potential is fixed by the scalar power 
spectrum Vs- In our model, from Eqs. ([20l) and ([T6| . it is clear that these observables depend on the transformation 
V pV. In this case, the transformations which keep a, (3, 7 and N constants are V pV and ^ ^/p. 



C. Reconstruction of potential and coupling functions 

Slow-roll inflation in the Gauss-Bonnet modified gravity violates the consistency relation r = —Sut- However, 
in the case of standard multi field inflation, this relation also modified as r < —Sut- Therefore, the case of over 
production of gravitational waves, r > — Srix, is important. From Eqs. (l68l) and (I69p . we can see over-production is 
realized if /3 > a is satisfied. Taking a look at the definition of slow-roll parameters and Eas. ((20l) and ((T4)) . we can 
see in the situation that > is realized, gravitational waves will always be overly produced. Thus, we can expect 
over production in many models. 

If the violation of this relation is observationally confirmed, we can reconstruct the value of the derivative of the 
Gauss-Bonnet coupling function From Eqs. ipO)) 

The potential V and the slow- roll parameter /3 are related by the scalar power spectrum Vs- From Ea. (H51) 

^^ = 24i^i4' ^^^^ 



Hence, ^_<i is written as 



Let us parameterize the violation ratio of consistency relation by using 6: 



S^'-±^. (74) 



The expression for becomes 



Mp.e. = 4^ = 1.86xl0^f Y—) • (75) 

^'^ TT^ y^Vs V2-441xlO^V VO.l/ VO-Oiy ^ ' 

This equation clearly shows that we can get the value of f^^, or constrain at least, with observing the tcnsor-to- 
scalar ratio r and the violation ratio of consistency relation 6. In other words, we can predict this violation ratio of 
consistency relation by using this equation for each Gauss-Bonnet infiationary model. The expression for the value 
of potential, at the CMB scale, is the same as usual. From Eq. ((72|) 



rVs = 3.61 X 10--^ 77T • (76) 



M^j 2 " ■ V2-441 X 10-9 J VO.l 

However, the derivative of potential suffers modification. From the definition of slow- roll parameters, Eq.(|2CT 

-10 f \ / r ^3/2 



4.04 X 10" 



2.441 X 10- 



/ r 



Using Eqs. ([75|) . ([7St and ([77]) . we can fix the value of V, and using observables. The expression for the scalar 
spectral index ng contains l^^^^ and ^,00. However, we cannot separate these quantities in our model. Thus, from 
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these observables, we cannot fix the value of Vcf,cf, and ^,0^. To this point, we have studied a general case, placing 
special emphasis on a violation of consistency relation. However, there are other models which also can violate this 
relation. Particularly, the model which overly produce tensor modes exists [s^l- Therefore in the next section, we will 
investigate an unique phenomenon in our model, namely a blue spectrum of tensor mode. 

Next, let us study the Chern-Simons coupling function. Considering future observations, deriving the expression 
for by using observable is meaningful. From the definition of il, Eq.(|52|). w,^ is 



uj 2M^p -1 _ 4M|i n _ 12Af]li n _ V2 n _ 16 V2 n 



--9^45X10M,,,;,^V m-'Y^I , (78) 



2.441x10-9/ VO.l/ VO-01 

Thus, we can see that for confirming existence of the Chern-Simons term, it is important to detect circular polarization 
in tensor modes. We can also use this expression to predict this polarization ratio. Because |n| <^ 1 is required by 
the definition. Eg. ([TO)) , we can get the constraint on the value of w^^: 



|Afp..,|«9.45xlO^(^^^^^^j . (79) 

Of course, if we allow an appearance of ghost in tensor modes, a situation might be changed drastically. However, 
it exceeds the scope of this paper. 



D. Lyth bound 

Here, we comment on the Lyth bound fss']. The statement of this bound is following: If primordial gravitational 
waves are detected, the variation of inflaton field must exceed a scale on the order of Planck mass. In our model, this 
bound also holds. From Eg. ([TO)) 

N=J-r^ = ^r^ (80) 
Mpi Jo VW Mpi Jo ^ ' 

and 

, Mpi Ir 



2 V2'^- ^''^ 

This can be easily understood: Large r corresponds to small 0, and this leads to small S(j). Note that effect of 
gravitational waves on the CMB multi pole is relevant in approximately I < 100, which corresponds to 6N ~ 4. This 
means that if gravitational waves with tensor-to-scalar ratio r are observed, the field value of inflaton has changed 
its value approximately \/2r during SN ~ 4. Hence, from this equation, we can set the lower bound on variation of 
inflaton field A0 during whole infiation: 

A0 ^ /— ^ / r \ i/2 



> V2r = 0.447 



PI 



/ r \ ^/-^ 

(al) ■ <«^' 



We can see that the Lyth bound holds in our case, while an equation of motion is modified due to the Gauss-Bonnet 
correction. Therefore, extremely small field inflation, namely Acj) <C Mp\, cannot produce an observable amount of 
tensor modes. 



IV. INTERESTING FEATURES OF OUR MODEL 



In our model, a spectrum of primordial gravitational waves can be blue or scale invariant, and these cases strongly 
violate the consistency relation r = — Snx. Hence, these (especially, a blue spectrum) are important for confirming 
our model with future observations. In this section, we use concrete examples to check these interesting situations 
are consistent with the current WMAP result. We show characteristics of potential and the Gauss-Bonnet coupling 
function which realize a blue or scale invariant spectrum. We also consider the case in which potential is cancelled 
out by the Gauss-Bonnet effective potential. 

We take e-folding number as iV = 60 for calculations in examples. 



11 



A. Blue spectrum 



The most interesting and unique phenomenon in our model might be a blue spectrum of tensor modes. Note that 
it is impossible to achieve a blue spectrum in a conventional model. For achieving a blue spectrum, a < is required. 
Hence < is necessary because 4> < 0. In a case of blue spectrum, a scalar field is climbing up potential, not 
rolling down. Of course, this blue situation must be over at some time, or inflation lasts forever, because H > is 
satisfied in this blue case. A climbing up situation can be understood using Ea. (IT51) . which shows that the coupling 
function ^ works as effective potential. Hence, if the effective potential ^ is steep enough, (j) can climbs up potential. 
We show this in FIGlH 

Vt 




FIG. 1: The configuration of the potential V and the coupling function ^, which achieve a blue spectrum, is shown. We need 
a climbing-up situation to realize a blue spectrum. Therefore the following situation is required: At early stage, the effective 
potential ^ makes (j> climb up potential and a blue spectrum is realized. At late stage, (j) rolls down potential as usual. 



In a blue spectrum case, can be calculated by employing Ea. (l75p with S > 1. Hence, there is a necessary 
condition for a blue spectrum: 



and combining with Ea. (|76p 



Remember S — 1 + Sn^/r and Inxl ^ 1- To observe tensor modes, r cannot take an extremely small value. Hence, 
5 cannot exceed a value of order unity and 5 1 must hold. In a consistent blue spectrum case, hence, we need the 
derivative of the Gauss-Bonnet coupling function on the order of 10^/Mpi and Mpj/F. It is important that both 
too small and too large a value of cannot result in a desirable blue spectrum. In other words, if a blue spectrum 
is detected, at the CMB scale must be on the order of 10*/Mpi. 

In the case that a blue spectrum of tensor modes is confirmed, we can reconstruct the value of ^,0, which must 
be on the order of 10^/Afpi. However, it is unclear whether or not an observationally consistent blue spectrum case 
can actually be realized; what kind of potential and the Gauss-Bonnet coupling function is likely to realize a blue 
spectrum. To show these, we use two types of examples of blue spectra. First one is a simple, but analytically solvable 
model and second one is a somewhat realistic model, which is calculated numerically. 

We take the following potential and coupling function: 



^^^{(j)) = const. (85) 
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In this case, the e-folding number N can be written analytically: 



,r2 2 

where 



71-2Af2j VI + 4612 i^x- ^2 - 1 



^=-^^™^ln(^^) , (86) 



TT 



Vbe,00o _ 1 + VI + 46*2 _ 1 - VI + 4612 ^ 1 



a; = cos , 6 = '■ — i— , xi = , xo — = . (87) 

\2<j>oJ ' SttM^j ' ^ 261 ' 29 xi ^ ' 

Note that x must be in the range — l<a;2<a;<l<a;i and 6 represents effect of the Gauss-Bonnet correction. We 
use a; as a variable instead of (j) and a; < corresponds to a climbing-up situation and a blue spectrum. To show a 
realization of blue spectrum, we express a; in ^ and N: 



1 + .T1 - (l + a;2)e2' VTTI^ + 1 + 261 + (VTTI^ - 1 - 26')e!' 

X ~ ^ — = I , , 88 

1 - X2 - {I - xi)ey y/TTW - 1 + 29 + (VTTW + 1 - 20)ey 

where 

y = "L^^^i + 46(2^ = y^^i + 45)2 (89) 
400 

Here, we have two parameters in this model, namely 9 and j/o- Let us consider varying to realize a blue spectrum 
in fixed yo, which corresponds to fixed (/(q. In the situation that 9^1, Eq. (|88| becomes 

45/ _ e^voS 

^ = A9 + e2yoe ' (9^) 

Hence, with sufficiently large 9, x < Q and blue spectrum can be realized for any (pQ value. The condition for a blue 
spectrum is 

A{l-q)9^ -Aqe+{l-q)<Q , (91) 

where 

g = tanh2(|) . (92) 

This inequality is satisfied, only if 1/2 < q < 1. We define t/b so that x = is realized &t y — y^- Note that y < j/b 
and y > yh correspond to red and blue spectra, respectively. Because 1/2 < g is necessary for a blue spectrum, the 
condition for j/b is 

2/b > yc = 2tanh~^ ^ = 1.76 . (93) 
v2 

Here, we define % and 6'b in the same manner as y^- From Eqs.dHH), (HH) and we can relate these variables: 



% = + 1 = tanh^ f ^) = tanh^ f , (g^^ 



(20b + 1)2 V 2 y V 

Note that for yo = yc/V2 model, 6'b — 1/2 (and yb = yc) is the solution for above equations. The scalar spectral 
index rig, the tensor spectral index rix and the tensor-to-scalar ratio r become 

ns-l=^-^--^[-2-x^ + 9x{l-x^)\ (95) 
iV 1 — a;^ 

n^ = ^^^[x + 9[l- X-)] = f Y^(x - -i)(- - -2) (96) 

' [x + 0(l-a:2)]^ = ^-ii^(x-.,)^(x-.2)^ (97) 



iV 1 - a;2 L ^ iV 1 - a;2 

Let us consider the large limit, in which y 3> yb is satisfied. In this limit, the following equation is realized 



2 



I+X1-X2+X2 , . 

x-X2 = -p. ^ . (98 

1 - X2 - (1 - a;i)eS' 
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Hence, x approaches to X2- We define new variable z: 

z EE X + e{l - x^) = -e{x - xi){x - X2) ■ (99) 

Because x ^ X2 is realized, z — > holds. Note that x — X2 approaches to zero exponentially as a function of 6, and 
this leads to 6{x — X2) — > 0. Observables become 

Note that r becomes small in large 9, because 9z — >■ 0. This shows that in any model, this large 9 limit results in 
small r and extremely red ng- The reason for small r is that x = X2 corresponds to the fixed point, at which potential 
is cancelled out by effective potential. Hence, = and r = are realized. Since this is not a desirable case, we 
consider the small 9 case, in which 9 ^ 9^ holds. In this case, a; ^ is realized, and observables become 

"s-l = -^, nx = — (-a;6ib) , = ^^^'b ■ (101) 

Note that 9\j and are not independent variables, because is fixed by choosing y^. Remember that 0b = 1/2 is 
the solution for ?/o = Vc/"/^ model. Let us substitute these: 

_ 1 = = -0.0415 , riT = ^ i-x) = -0.0104X , r = = 0.0415 . (102) 

N 2V2N ' N 

Here, we get good result: ns and r are consistent with current observations and a blue spectrum of tensor modes is 
obtained. Therefore, desirable solutions might exist in the vicinity of ^ 1/2 and yo ^ yc/%/2: 

0o-^^Mp> = lO.9Afp,, ^-^/l- 0.432. (103) 

Let us consider the model, in which (po is large enough to satisfy yo ^ Vc- From Eas. ([551) and yh — '2yo9h is 

realized, and r at a; = becomes 

r = ^^byb > ^^b . (104) 

Because 6'b = yb/2yo > J/c/2?/o ^ 1, this large field model results in r, which is too large. This is because, in this large 
field model, infiaton (j) must roll large field value between the hill-top of potential and the vacuum, within iV = 60 
e-fold. This requires large 0, in other words, large r. In addition, from Eg. dlOip . this large field model also leads to 
a flat spectrum of scalar modes. In the small field model, which satisfies yo S> yc, the inequality y ^ 1 always holds. 
Hence, as in the above large 9 limit, this results in small r and extremely red ns. Note that larger field models are 
corresponding to larger values of r, and this is consistent with the consideration in the Lyth bound. 

We show 6'-trajectories of solutions on a ns — r plane in FIGlH We plot the 68% and 95% C.L. contours from 
the WMAP 7-year result[l|. We choose eight (po values surrounding 00 = llA/pi. Thick dotted and thick solid 
lines correspond to 0.003 > nx > and nx > 0.003 blue spectra. As stated above, models with 0o ^ llMpi 
have observationally consistent blue regions; large field models correspond to large r and small or no consistent blue 
spectrum region; small field models have extremely small r; in large 9 limit, r becomes zero and ns becomes extremely 
red. It is clear that a consistent blue spectrum model can actually be constructed. To get the value of ^^0, we can use 
Eg. ([75)1 . Of course, as mentioned above, ~ 10^/Mpi is realized in a blue region. 

Next, we consider a realistic model. To realize a blue spectrum, new-infiation-type potential, as depicted in FIG[1] 
is required. Therefore, we use forth-order double-well potential and an exponential coupling function: 



V{<j)) = plO-i°Mj*i 



^U) ^ s- exp(sa0/Mpi) , (105) 
P 



where s — ±1, p will be fixed by the scalar power spectrum Vs = 2.441 x 10~^[l| and a,b > 0. We take (p — as our 
vacuum. We use this exponential coupling, since it appears in many theories, such as a dilatonic case or an one-loop 
correction from a superstring motivated model[2]. 

We show 6-trajectories on a ns — r plane in FIG|31 We plot the 68% and 95% C.L. contours from the WMAP 7- year 
result. We use nine values of a, which varies from 0.01 to 1 with an even logarithmic interval, and three value of (f>o, 
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ns 

FIG. 2: The trajectories, in which d is varying, of analytic blue models are shown on a ns — r plane. Two contours roughly 
corresponds to the 68% and 95% confidence level of the WMAP 7-year result. We draw eight lines, which are corresponding 
to <^o = 9 (red, bottom) to 16Mpi (blue, top) with an even interval. Thick dotted and thick solid lines denote the blue spectra 
of 0.003 > riT > and tit > 0.003, respectively. Models with 00 ~ llMpi are consistent with observations and realizes blue 
spectra. 



11.1, 15 and 20Mpi. We choose 0o value so that ns and r values are inside this 68% C.L. contour in the 6 — > Umit. 
Thick dotted and thick sohd hues in FIGl3] denote the blue spectra of 0.003 > nx > and nx > 0.003, respectively. 
We can see it is possible to realize a consistent blue spectrum. The case of 0o — ISAfpi is appropriate to achieve a 
consistent blue spectrum. Small (^^ leads to small r and large 0o leads to a small consistent blue region, and these are 
essentially consistent with the above analytic case. Due to this exponential coupling, even large field models, namely 
(^0 = 20Mpi models with s = —1, have blue spectrum regions and small field models, namely (pQ = ll.lA/pi models 
with s = +1, have an observable amount of r. This figure shows that an order unity value of a leads to a non-desirable 
result, which has too small r or no blue region. Therefore a <C 1 is appropriate to obtain a blue spectrum. Moreover, 
if a is large in s = —1 models, it implies the Gauss-Bonnet effect becomes stronger toward the end of inflation. In 
this case, standard reheating scenarios seem to be broken. Therefore, large a in s = — 1 models is also not desirable 
in this sense. Note that with Ea. ([75)) . we can get the value of 

From the above analytic model and this realistic one, we can say that inflationary models with observationally 
consistent blue spectra can be constructed. Hence, the Gauss-Bonnet coupling function will actually be fixed or 
constrained with future observations. New-inflation-type potential with the symmetry breaking scale of (j)o ~ lOMpi 
and an almost constant form of might lead to a consistent blue spectrum. Of course, these are not general 
conditions. However, we can expect many models are also consistent with these conditions. Note that to obtain blue 
spectrum, ~ 10*/A/pi must be realized in any model at the CMB scale. 



B. Scale invariant spectrum 



In inflation with a standard gravitational action, a scale invariant spectrum of gravitational waves corresponds to 
almost purely de Sitter inflation and this requires the tensor-to-scalar ratio r extremely small. However, inflation with 
the Gauss-Bonnet term behaves differently. Let us consider the case, in which the Gauss-Bonnet term dominate the 
scalar equation of motion, Eq. p5|) . The function ( becomes 



C = , (106) 

^ 6M4 ' ^ ' 



PI 



and this requires 



^ = - » 1 • (107) 



In this case, a becomes negligible. Hence 

ns - 1 = 27 = 2M^M , nx = , r = 16^ = 8M|X . (108) 
V V 
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FIG. 3: The trajectories, in which b is varying, of some blue models are shown on a ns — r plane. In the 6^0 limit, all 
trajectories converges to one point. Two contours roughly corresponds to the 68% and 95% confidence level of the WMAP 
7-year result. We draw nine lines, which are corresponding to a = 0.01 (red, left in left panels and top in right panels) to 1 
(blue, right in left panels and bottom in right panels) with an even logarithmic interval. Thick dotted and thick solid lines 
denote 0.003 > wt > and ut > 0.003 blue spectra, respectively. The model with (po ~ 15A/pi is appropriate to achieve blue 
spectrum. An exponential coupling makes large (po i4'o = 20A/pi) model result in a blue spectrum, and small (po {4>o = ll.lMpi) 
model achieve an observable amount of r. 
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These equations means that if the Gauss-Bonnet correction dominates this equation of motion, the tensor-to-scalar 
ratio r can take an observable value although rix = is realized. Because of the consistency relation r = — Snx, this 
scale invariant spectrum is impossible in a conventional model. Therefore, if scale invariant primordial gravitational 
waves are detected, it implies existence of the Gauss-Bonnet term in an inflationary theory and is determined by 
Eq.dZSl) with 5=1: 



and using Eg. ([75)) 



V2 1 



1.86 X 10*^ 



Vs 



2.441 X 10- 



r 

OA 



-1/2 



1/2 



(109) 



(110) 



As in a blue case, ~ 10^/Mpi and ~ Mp^/V must be satisfied in this flat spectrum case, too. Both too large 
and too small a value of cannot lead to a consistent scale invariant spectrum. From Eas. (|107p and (|108p . the 
inequality for r holds: 



vi 



(111) 



This shows that potential must be nearly flat for the tensor-to-scalar ratio r to remain small. 

Here, we check whether or not an observationally consistent scale invariant model can be constructed. In this 
subsection, we also use two types of examples, namely a simple and analytic one and a realistic one. Let us start with 
an analytic model: 



V = Vo 



tanh 



, C,0(0) = const. 



(112) 



The e-folding number TV becomes 



N = 



2M2j 



(113) 



where x and are defined by 



X = sinh I — 



^ VoL<i>'f>o 



(114) 



We use X instead of (j), and x is 



exp 



(115) 



The observables, namely the scalar spectral index ng, the tensor spectral index and the tensor-to-scalar ratio r 
become 



ns - 1 = 



riT 



Ml, 1 ( 


(t>l x"^ + 




Ml, 1 


t(- 


X? + 




%Ml, 1 









X^ 

1 

x^ 



6x 



(116) 
(117) 
(118) 



In this model, we also variate 9 and fix 

Let us consider the limit case, in which 9 <^ (/)g/2A/pjiV is satisfied, and x becomes 



(119) 
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Hence, observables are 

1 4a;2 + 3 11 4 1 , ^ 

We can see that in this small 9 limit, all models give observationally consistent values. Next, we investigate the large 
6 limit, ^ 0Q/2AfpjiV. The variable x takes an exponentially large value in this limit. Hence, observables become 

These show that in this limit, ng and nx approach to flat one and r becomes large. Remember that x behaviors 
exponentially as a function of 9. If this large 9 limit is realized in the WMAP consistent region, such as r < 0.3, the 
result which has an observable amount of r and flat nx is achieved. Here, we consider 9 — (ty^/2M'^^N case, which is 
the boundary case between large and small 9. In this case, 1 + x'^9 = e is realized. The observables nx and r become 

I e 9 9 A 9 9 , , 

nx = = -0.0132 , r = = 0.287 . (122 

2Afe-le-l + 6i c-l + 6i iVc-le-1 + 6' 6-1 + 6* ^ ' 

li 9 = (j)l/2MpiN > e - 1 is satisfied, -nx > 0.00659 and r > 0.143 hold. This r might be too large to reahze that 
large 9 limit in r < 0.3. Therefore, to achieve desirable results, 6* ^ c — 1 at this 9 = <pQ/2MpyN case might be 
necessary. In other words. 



(/>o < y'2N{e - l)Mpi = 14.4Afpi . (123) 

The physical meaning of this condition can be easily understood: Large field model requires large (j) and inconsistently 
large r(> 0.3), to reach a flat potential region within = 60 e-fold. This interpretation is essentially the same as 
that in the condition for blue spectrum. The difference is that small field model can also achieve large r in this case, 
because small doesn't means the smallness of the variation of the field value of the infiaton (j). It means the smallness 
of the field value between the vacuum and a flat region of potential. Hence, this is consistent with the Lyth bound. 
If this condition is satisfied, sufficiently large 9, namely 9 3> (j)Q/2Mp^N , in other words 

Ma » (124) 

M|i ^ N Mpi ' ^ > 

leads to a consistent fiat spectrum. 

We show ^-trajectories of solutions on a ns — r plane in FIG|H We plot the 68% and 95% C.L. contours from 
the WMAP 7-year result. We choose five (po values, 2, 4, 6, 8 and lOMpi. Thick lines denote the fiat region, 
in which r/8|nx| > 100 is realized. Note that r = — 8nx holds in conventional single field inflation. Therefore 
large r/8|nx| value implies a flat spectrum. We can see that small (po models, namely ipo — 2Mp\ and ipo = 4A/pi 
models, result in WMAP-consistent flat spectra. This is consistent with the criterion at the above calculation, namely 
(j)o ^ 14.4Afpi. As mentioned above, ns also becomes flat, in a flat tensor region. Note that the smaller fleld model, 
namely 0o — 2A/pi model, has larger consistent scale invariant region. Of course, ^.^ can be reconstructed from 
Ea. (l75|) . and ^.^ ^ lO^/A/pi must hold in a flat region. 

Here, we consider a realistic example of a flat spectrum model. Nearly flat potential results in a Gauss-Bonnet 
dominating situation. Therefore we take the following potential and coupling function: 



V{(l)) = plO^i^Af^i 



1 — exp — -pj 



/,2 



^(0) = s- exp(sa0/Afpi) , (125) 
P 



where s = ±1, p will be flxed by Vs — 2.441 x 10^^, and a,b > 0. In the — oo limit, this potential becomes 
extremely flat. We use the same exponential coupling in this example as in previous subsection. 

We show ^-trajectories on a ns — r plane in FIGl5] We plot the WMAP 68% and 95% C.L. contours. We use 
three values of (/)o, 2, 6 and lOA^pi, and the same set of a in this case as in the realistic blue spectrum case. Thick 
lines denote the regions r/8|nx| > 100. As predicted from the previous analytic models, smaller 4>o models has larger 
regions of flat spectra. Note that s — —1 models have larger consistent flat region than that of s = — 1 models. An 
order unity value of a leads to extremely small r in both s — ±1 models, and considering reheating, large a in s = — 1 
models is not desirable. Therefore also in this case, a ^ 1 is appropriate. Eq. ([75|) is valid to get the value of ^,0. 
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FIG. 4: The trajectories, in which 9 is varying, of analytic flat models are shown on a ns — r plane. Two contours roughly 
corresponds to the 68% and 95% confidence level of the WMAP 7-year result. We draw five lines, which are corresponding to 
(j>o = 2 (red, bottom) to lOMpi (blue, top) with an even interval. Thick lines denote the regions, which realize r/|8nT| > 100. 
Small (po models, namely (j>o ~ 2Mpi and (j>o = 4A/pi models, realize observationally consistent flat spectra. 



In s = —1 models, trajectories approach to some points for each a in large b limit, in FIGE] This is understood 
from Eg. dTSl) . In the large b limit of s = — 1 models, only a large 4> region contribute to the integration for an e- folding 
number. Because V<p/V ~ in this region, we can omit V<p/V term: 

iV=^r^.4^. (126) 



Vo Jo aVo^,, 

This shows that, in this large b limit, approaches to ^.^ = 6A/p[/aA^Vo. This also fixes the value of ^.^^ and all 
slow-roll parameters. The observables ns and r become 

, 2 2Vo^e^ 8 
ns — 1 = -TT- = , r — TT-^ = — — - . (127 

Hence, trajectories approach to some points in large b. This shows that a ~ 1 model leads to very small r because r 
approaches to r = = 0.00222 at a = 1. 

From above results, we can say that an observationally consistent flat spectrum can be obtained. Hence, confirming 
existence of the Gauss-Bonnet term is possible, with observing a scale invariant tensor mode. In the cases we 
considered, (j)o <^ lOA/pi and nearly constant ^,0 might lead to a consistent flat spectrum. Note that (j)o denotes the 
difference between the filed value of a flat region and the vacuum of potential. We can expect these conditions are 
also satisfied in many other cases. To realize a flat spectrum, 10^/Afpi is necessary in any model. 



C. Inflation with steep potential 



As mentioned in the end of section |TT1 if we take ~ 6Mp^/V, potential is cancelled out by the Gauss-Bonnet 
effective potential. Considering this fact, let us take ^ as follows: 

f^(l-»)5f^. (128) 

where k is a constant and we assume \k\ <C 1. The function becomes 

C = - (1 - = • (129) 

Eg. dTSI) becomes 
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FIG. 5: The trajectories, in which b is varying, of some flat models is shown on a ns — r plane. In the 6^0 limit, all trajectories 
converges to one point. Two contours roughly corresponds to the 68% and 95% confidence level of the WMAP 7-year result. 
We draw nine lines, which are corresponding to a = 0.01 (red, left in left panels and right in right panels) to 1 (blue, right in 
left panels and left in right panels) with an even logarithmic interval. In this case, tensor spectrum is always red. Thick lines 
denote the regions of flat spectra, r/8|nT| > 100. This figure shows that small 00 leads to a fiat spectrum. Note that in s = — 1 
models, trajectories approach to some points, in large b limit. 



20 



Hence if k < is realized, a scalar field always rolls up potential toward a hill-top. This is not a desirable situation. 
Therefore we take < k ^ 1. Slow- roll parameters becomes 

P = -'^^^0, 7 = «M|,^, (131) 

where we used k <C 1. In this case, /3 becomes negligible. Note if the Gauss-Bonnet term is dominating, a becomes 
negligible. In this case, the observables n^, and r are 

ns - 1 = -6a + 27 , ht ^ -2a , r = 0. (132) 

Because the tensor-to-scalar ratio r is zero, the value of the spectral index ut is meaningless. Therefore observationally 
speaking, this situation is not interesting. However, because k 1, three slow- roll conditions are automatically 
satisfied. Hence, even if we use extremely steep potential, inflation can be achieved. 



V. CONCLUSION 



In this paper, we studied slow-roll inflation with the Gauss-Bonnet and Chern-Simons corrections. We defined 
three slow-roll parameters, while we used five parameters in the previous paper 35]. We derived expressions for the 
scalar spectral index ng, the tensor spectral index the tensor-to-scalar ratio r and the circular polarization ratio 
H by using these parameters. We showed that in our model, the consistency relation r = —8n^ is not automatically 
satisfied. If this violation is observationally confirmed, we can determine the derivative of the Gauss-Bonnet coupling 
function at the CMB scale. In our model, even both blue and scale invariant spectra of gravitational waves can 
be realized. Because blue and scale invariant mean nx > and |8nT|/f ^ 1, respectively, these cases violate this 
consistency relation strongly. Therefore these are the key for confirming our model in future observations. We showed 
that if either a blue spectrum or a scale invariant one is observed, it supports existence of the Gauss-Bonnet coupling 
function of ^ ^ 10^/Mpi at the CMB scale. We checked whether or not these blue and scale invariant spectra are 
consistent with current observations. For this purpose, we used concrete examples. We showed that new-inflation- 
type potential with lOMpi symmetry breaking scale and potential with flat region in > lOMpi might result in 
observationally consistent blue and scale invariant spectra, respectively. An almost linear form of the Gauss-Bonnet 
coupling function is appropriate in both cases. These tell us that the detection of a blue or scale invariant spectrum 
of tensor modes can be actually expected with future observations. We also showed that if circular polarization of 
gravitational waves is detected, the derivative of the Chern-Simons coupling must be on the order of 10*/Mpi. 
Thus, we can say that existence of gravitational higher coupling terms in a inflationary model will be confirmed, or 
at least constrained, with future experiments. 

These higher curvature terms are expected to appear in non-Gaussian part of perturbations. Therefore, calculating 
non-Gaussianity in our model is important future work. 



We would like to thank Jiro Soda for his very useful comments. 
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